Spin Dynamics of the Spin-1/2 Kagome Lattice Antiferromagnet ZnCu 3 (OH) 6 Cl2 
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We have performed thermodynamic and neutron scattering measurements on the S — 1/2 kagome 
lattice antiferromagnet ZnCu3(OH)6Cl2. The susceptibility indicates a Curie- Weiss temperature of 
6c w — —300 K; however, no magnetic order is observed down to 50 mK. Inelastic neutron scattering 
reveals a spectrum of low energy spin excitations with no observable gap down to 0.1 meV. The 
specific heat at low-T follows a power law temperature dependence. These results suggest that an 
unusual spin liquid state with essentially gapless excitations is realized in this kagome lattice system. 
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An important challenge in condensed matter physics is 
the search for quantum disordered ground states in two 
dimensional systems. Of particular interest is studying 
quantum spin liquids, an example of which is the "res- 
onating valence bond" state proposed by Anderson [lj]. 
These states are unusual in that neither translational nor 
spin rotational symmetries are broken. It is believed that 
the 5=1/2 Heisenberg antiferromagnet on a kagome lat- 
tice (composed of corner sharing triangles) is an ideal sys- 
tem to look for spin liquid physics due to the high degree 
of frustration. There is broad theoretical consensus that 
the ground state of the 5 = 1/2 kagome antiferromagnet 
is not magnetically ordered [3, d, [1, @, 0, [!]■ How- 
ever, many basic properties are still under debate, such 
as the magnitude of the gap to the first triplet state. 
An intriguing possibility is the existence of deconfmcd 
5 = 1/2 spinons as the fundamental excitations, as op- 
posed to conventional 5=1 magnons. 

Despite heavy theoretical interest, experimental stud- 
ies of the 5 = 1/2 kagome lattice have been ham- 
pered by the difficulty in synthesizing such materials. 
Here, we report thermodynamic and neutron scattering 
measurements on powder samples of ZnCu3(OH) 6 Cl2, 
known as herbertsmithitej^]. As has been previously 
reported[I3|, Zn a; Cu4_ a; (OH) 6 Cl2 can be synthesized 
with variable Zn concentration, from x—0 to x—1 (her- 
bertsmithite). Figure 1(a) represents the transfor- 
mation from Cu2(OH)3Cl, which has a distorted py- 
rochlore structure, to ZnCu3(OH) 6 Cl2, which consists of 
Cu kagome layers separated by nonmagnetic Zn layers. 
Structurally, ZnCu3(OH)gCl2, with space group i?3m 
and lattice parameters a = b = 6.832 A and c = 14.049 A, 
appears to be an excellent realization of the 5 = 1/2 
kagome lattice antiferromagnet. Initial evidence is the 
absence of long-range magnetic order, as shown in the 
neutron diffraction scans in Fig. 1(b). In Cu2(OH)3Cl, 



clear magnetic Bragg peaks are observed below ~ 6 K; 
whereas no magnetic Bragg peaks are observable down 
to 1.8 K in ZnCu 3 (OH) 6 Cl 2 . 

To further characterize the properties of 
ZnCu3(OH)gCl2, we performed magnetic susceptibility 
measurements on powder samples. The susceptibility, 
shown in Fig. 1(c), can be fit to a Curie- Weiss law at high 
temperatures (T > 200 K). The resulting Curie- Weiss 
temperature of —300 ±20 K implies an antiferromagnetic 
exchange J ~ 17 meV, calculated using the series expan- 
sion corrections for the kagome lattice 11, 12, fl3| . The 



susceptibility continually increases as the temperature 
is lowered down to 1.8 K. At first glance, this behavior 
may suggest the presence of several percent free spin-1/2 
impurities yielding a Curie tail. This is certainly possi- 
ble, but is not necessarily the case. From the chemical 
analyses, we calculate the stoichiometric coefficients to 
be 3.00 ± 0.04 on the Cu site and 1.00 ± 0.04 on the Zn 
site. Also, we have measured the ac susceptibility at 
temperatures down to 50 mK, as shown in the inset of 
Fig. 1(c). These data do not follow the simple Brillouin 
function behavior expected for free 5 = 1/2 spins. In 
particular, the susceptibility increase from 705 mK to 
50 mK is much smaller than the free spin prediction. 
Recently, Ofer and coworkers [l^ have shown that the 
muon Knight shift and transverse relaxation rate have 
T dependences similar to the measured susceptibility. 
Hence, the measured susceptibility may be intrinsic to 
the Cu kagome system. We note that similar behavior 
is found for the frustrated 5 = 1/2 nuclear moments of 
3 He films on graphite, where the susceptibility is found 
to continually increase with decreasing temperature 
down to T ~ J/300 [H. Another recent /xSR study [H 
emphasizes the role of defects. The roles of impurities 
and exchange or Dzyaloshinskii-Moriyafl^ anisotropies 
in this system remain important topics for further 
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FIG. 1: (color online) (a) The chemical transformation from 
the pyrochlorelike lattice of Cu2(OH)3Cl to the kagome lay- 
ers of ZnCu3(OH)6Ci2. (b) Magnetic diffraction scans of the 
two systems at T — 1.4 K (open) and 20 K (filled). The 
Cu2(OH)3Cl data show magnetic Bragg peaks at Q ~ 0.70 
and Q ~ 0.92 which are absent for the ZnCu3(OH)6Ci2 data 
(which have been shifted by 2300 cts/min for clarity), (c) 
Magnetic susceptibility of ZnCu3(OH)6Ci2 measured using a 
SQUID magnetometer plotted as 1/x, where mole refers to a 
formula unit. The line denotes a Curie- Weiss fit. Inset: ac 
susceptibility (at 654 Hz) at low temperatures measured at 
the NHMFL in Tallahassee, FL. 

investigation. We also observe a small peak in the ac 
susceptibility near H — 2 T at 50 mK which disappears 
upon warming to 705 mK. The overall susceptibility 
data indicate the absence of magnetic order or a spin 
gap down to 50 mK. 

The specific heat C(T) of ZnCu 3 (OH) 6 Cl 2 is shown in 
Fig. 2(a) in various applied fields. For temperatures of 
a few Kelvin and higher, the lattice contribution to the 
specific heat (proportional to ~ T 3 ) is the most signif- 
icant contribution, as shown in the inset. However this 
contribution diminishes at low temperatures, and below 
~ 5 K, an additional contribution is clearly observed 
which arises from the Cu spin system. Magnetic fields 
of a few Tesla can significantly affect the low-T behavior, 
and fields of 10 Tesla and higher strongly suppress the 
specific heat below 3 K. The difficulty in synthesizing an 
isostructural nonmagnetic compound makes it hard to 
subtract the lattice contribution precisely. However, the 
magnetic field dependence suggests that the specific heat 
in zero applied field below 1 K is predominately magnetic 
in origin. As a rough measure of the spin entropy, the 
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FIG. 2: (color online) (a) The specific heat C(T) of 
ZnCu3(OH)eCl2 in various applied fields, measured using a 
Physical Properties Measurement System. Inset: C(T) plot- 
ted over a wider temperature range in applied fields of T 
(square) and 14 T (star), (b) C(T) in zero field measured 
down to 106 mK. The lines represent power law fits as de- 
scribed in the text. 

field-induced change in specific heat below 3 K, obtained 
by subtracting the 14 T data from the zero field data, 
accounts for about 5% of the total entropy of the spin 
system. 

Additional specific heat measurements at zero field at 
temperatures down to 106 mK were performed at the 
National High Magnetic Field Laboratory (NHMFL) and 
the combined data are shown in Fig. 2(b). The specific 
heat at low temperatures (T < 1 K) appears to be gov- 
erned by a power law with an exponent which is less than 
or equal to 1. In a 2D ordered magnet, magnon excita- 
tions would give C ~ T 2 . The kagome-like compound 
SrCr 8 _ a; Ga4+ a; Oi9 (SCG0)[3 and other 2D frustrated 



magnets [19[ are also observed to have C 
absence of long-range order [2< 



T even in the 



21| . The behavior that we 



observe in ZnCu3(OH) 6 Cl2 below 1 K stands in marked 
contrast. We can fit our data to the power law C — -fT a , 
though we note that the exponent a is sensitive to the 
chosen range of temperatures that are fit. The blue line in 
this figure represents a linear fit with a = 1 over the tem- 
perature range 106 mK < T < 400 mK. The fitted value 
for 7 is 240 ± 20 mJ/K 2 Cu mole. If we include higher 
temperatures, the red line represents a fit with a = 2/3 
over the temperature range 106 mK < T < 600 mK. Ex- 
tending the fitted range to even higher temperatures can 
yield a values as low as 0.5. 

Finally, inelastic neutron scattering measurements of 
the low energy spin excitations were performed on deuter- 
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ated powder samples of ZnCu3(OD) 6 Cl2. High resolu- 
tion measurements were taken on the time-of-flight Disk 
Chopper Spectrometer (DCS) at the NIST Center for 
Neutron Research in Gaithersburg, MD. A sample with 
mass 9 g was cooled in a dilution refrigerator and stud- 
ied with incident neutrons of wavelength 7 A, yielding an 
instrumental energy resolution of 0.02 meV (half-width 
at half-maximum). As shown in Fig. 3(a), the spin ex- 
citations form a broad spectrum at low energies. A no- 
table observation is the near temperature independence 
of the scattering for positive energy transfers. The ex- 
citation spectrum on the negative energy-transfer side is 
suppressed at low temperatures due to detailed balance. 

The magnetic scattering intensity is proportional 
to the dynamic structure factor S(Q,u>) = {n{uj) + 
l)x"(Qi w )j where n(u) is the Bose occupation factor and 
x"(Q, w) is the imaginary part of the dynamical suscep- 
tibility. We find that part of the measured intensity for 
positive energy transfers below 0.4 meV is spurious back- 
ground scattering, probably caused by multiple scatter- 
ing of neutrons within the sample environment. There- 
fore, to extract the intrinsic scattering from the sample, 
the following procedure was used. For negative energy 
transfers, x"(w,T= 10K) can be obtained by subtract- 
ing the 35 mK data (which is essentially background) 
from the 10 K data and dividing by the Bose factor. 
Here, x"( w ) represents the dynamical susceptibility in- 
tegrated over momentum transfers 0.25 < \Q\ < 1.5 A -1 
and is a good measure of the local response function. 
This is plotted in Fig. 3(b), where the positive ui data 
is obtained by using the fact that x"(w) is an odd func- 
tion of lo. Then, x"(cj,T= 35mK) can be extracted from 
the positive energy transfer data using S(u>; T— 35mK) — 
S(u; T= 10K) = I(uj; T= 35mK) - I(uj; T= 10K), where 
I(ui) is the measured intensity and the background is as- 
sumed to be temperature independent between 35 mK 
and 10 K. As seen in Fig. 3(b), the data for at 
T = 35 mK increase with decreasing w, indicating the 
absence of a spin gap down to 0.1 meV. Moreover, the 
data may be described by a simple power law; the solid 
line represents a fit to the form x"( w ) °^ tj7 with an expo- 
nent 7 = —0.7 ±0.3. This apparently divergent behavior 
is unusual and again differs markedly from measurement 
on SCGOji^] which yield 7 ~ 0. Of course, within the er- 
rors, we cannot rule out other functional forms for x"(ui). 

The Q-dcpcndence of the scattering is shown in 
Fig. 3(c). These data were obtained by integrating over 
energy transfers —0.5 < hcu < —0.22 meV and subtract- 
ing the 35 mK data set from the 10 K data set. We find 
that the data appear to be only weakly dependent on \Q\. 
Note that due to the polycrystalline form of the sample, 
the data represents the powder average of the scatter- 
ing from a crystal. The solid line represents the squared 
form factor |F| 2 for the Cu 2+ ion. The deviations of the 
data from \F\ 2 suggest that the structure factor is not 
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FIG. 3: (color online) (a) Inelastic neutron scattering data 
taken on DCS, integrated over momentum transfers 0.25 < 
|Q| < 1.5 A -1 , (b) x'{uj), extracted from the data as de- 
scribed in the text. The line denotes a power law fit. (c) 
The Q-dependence of the scattering, integrated over energy 
transfers —0.5 < tuv < —0.22 meV. (d) Energy scans taken 
on SPINS at zero field and 11.5 T at \Q\ = 0.6 A" 1 and 
T — 1.2 K. The lines are guides to the eye. Inset: Temper- 
ature dependence of the scattering for 0.3 < %u) < 0.5 meV 
and |Q| = 0.9 A" 1 . The blue data point and line indicate the 
background, measured on the energy loss side at T = 1.5 K. 

completely independent of \Q\. That is, some degree of 
spin correlations are necessary to account for the relative 
reduction in scattering at small \Q\. The overall diffuse 
nature of the scattering points to the absence of a well- 
defined length scale to describe these correlations. 

Further measurements were taken using the triple- 
axis SPINS spectrometer at the NCNR with the sample 
mounted inside a superconducting magnet, as shown in 
Fig. 3(d). The instrument was configured in a horizon- 
tally focusing analyzer geometry with Ef — 3.05 meV 
and collimations of guide-80'-radial-open. A BcO filter 
was placed in the scattered beam to reduce higher-order 
neutron contamination. The resulting instrumental en- 
ergy resolution was about 0.06 meV. An applied field of 
11.5 Tesla transfers spectral weight from lower to higher 
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energies. This demonstrates that a significant fraction of 
the low energy scattering is magnetic in origin, since the 
incoherent and phonon background would not respond 
to the applied field in this manner. The magnetic signal 
in zero field extends down to below 0.2 meV, consistent 
with the analysis of the above DCS data. In 11.5 T, the 
magnetic signal becomes peaked around Tito ~ 1.4 meV, 
which is close to the Zeeman energy g^sH. However, the 
half-width of this peak of about 0.21 meV is significantly 
broader than the resolution. Therefore, the peak does not 
simply originate from Zeeman excitations of noninteract- 
ing spins, which would result in a narrow energy peak, 
but involve spins which are part of the interacting system. 
The integrated spectral weight of the zero field magnetic 
signal for Kw < 1 meV accounts for at most 20% of the 
total scattering expected from a, S — 1/2 spin system (an 
estimate made by normalization to the incoherent elas- 
tic scattering from the sample and also to a vanadium 
standard). The inset of Fig. 3(d) shows the temperature 
dependence of the inelastic signal with energy transfers 
integrated over the range 0.3 <Tiu>< 0.5 meV. There is 
a small increase in the signal when the sample is cooled 
below ~ 5 K, though, for the most part, the intensity is 
largely independent of temperature in this range. 

Our experimental results suggest an intriguing pic- 
ture for the ground state properties of the 5 = 1/2 
kagome lattice antiferromagnet. A hallmark of the quan- 
tum spin liquid in 2D is the existence of deconfincd 
5=1/2 spinons as the fundamental magnetic excita- 
tion. A rich variety of spin liquid states have been theo- 
retically proposed in which the spinons can be described 
as bosonicja H, EH , fermionic 24 , [2|| , or even as Dirac 
fermions[26|. We note that several of these theories are 
based on triangular lattice Hamiltonians, and they may 
not have clear extensions to the kagome lattice antiferro- 
magnet. Using a naive comparison to a generic model of 
fermionic spinons with a Fermi surface, one would expect 
C = jT. From our linear fit below 400 mK, the value 
of 7 indicates a Fermi temperature of Tp ~ 110 K. How- 
ever, other forms for the specific heat (such as C oc T 2 ) 
may hold at higher temperatures where the lattice contri- 
bution prevents us from clearly identifying the magnetic 
contribution. 

The neutron scattering measurements of the excitation 
spectrum at low temperatures are also consistent with ex- 
pectations of deconfined spinons in a spin liquid. We find 
no evidence of a spin gap down to ~ J / 170, much lower 
than the prediction from exact diagonalization studies for 
a spin gap of - J/200. The power law behavior of x"( w ) 
is interesting and may indicate a spin liquid with criti- 
cal spin correlations [27J. Our observation of a diffuse Q 
dependence for the inelastic scattering suggests that if a 
singlet spin liquid picture is correct, then the singlets are 
not restricted to nearest neighbor dimers, since no well- 



defined length scale is indicated by the data. The near 
temperature independence of S(Q,lo), similar to obser- 
vations in /-electron systems (Hf, may indicate the prox- 
imity to a quantum critical point. Many of the current 
theories for 2D spin liquids were formulated to describe 
experimental results [2!| H(| for 5=1/2 triangular lat- 
tice systems. More theoretical studies based explicitly 
on the 5=1/2 kagome Heisenberg antiferromagnet (in- 
cluding the possible effects of impurities and exchange or 
Dzyaloshinskii-Moriya anisotropies) are certainly impor- 
tant for further comparisons with experimental results. 
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